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ABSTRACT 

We investigate the viscous two temperature accretion disc flows around rotating black holes. 
We describe the global solution of accretion flows with a sub-Keplerian angular momentum 
profile, by solving the underlying conservation equations including explicit cooling processes 
selfconsistently. Bremsstrahlung, synchrotron and inverse Comptonization of soft photons 
are considered as possible cooling mechanisms. We focus on the set of solutions for sub- 
Eddington, Eddington and super-Eddington mass accretion rates around Schwarzschild and 
Kerr black holes with a Kerr parameter 0.998 . It is found that the flow, during its infall from the 
Keplerian to sub-Keplerian transition region to the black hole event horizon, passes through 
various phases of advection - general advective paradigm to radiatively inefficient phase and 
vice versa. Hence the flow governs much lower electron temperature ~ 10 8 - 10 9 5 K, in the 
range of accretion rate in Eddington units 0.01;M;100, compared to the hot protons of tem- 
perature ~ 10 10 2 - 10 11 8 K. Therefore, the solution may potentially explain the hard X-rays 
and y-rays emitted from AGNs and X-ray binaries. We then compare the solutions for two dif- 
ferent regimes of viscosity and conclude that a weakly viscous flow is expected to be cooling 
dominated, particularly at the inner region of the disc, compared to its highly viscous counter 
part which is radiatively inefficient. With all the solutions in hand, we finally reproduce the ob- 
served luminosities of the under-fed AGNs and quasars (e.g. Sgr A*) to ultra-luminous X-ray 
sources (e.g. SS433), at different combinations of input parameters such as mass accretion 
rate, ratio of specific heats. The set of solutions also predicts appropriately the luminosity 
observed in the highly luminous AGNs and ultra-luminous quasars (e.g. PKS 0743-67). 

Key words: accretion, accretion disc — black hole physics — hydrodynamics — radiative 
transfer 



1 INTRODUCTION 

The cool Keplerian accretion disc (Pringle & Rees 1972; Shakura & Sunyaev 1973; Novikov & Thorne 1973) was found to be inappropriate to 
explain observed hard X-rays, e.g. from Cyg X-l (Lightman & Shapiro 1975). It was argued that secular instability of the cool disc swells the 
optically thick, radiation dominated region to a hot, optically thin, gas dominated region resulting in hard component of spectrum ~ lOOKeV 
(Thorne & Price 1975; Shapiro, Lightman & Eardley 1976). This region is strictly of two temperatures with electron and ion temperatures 
respectively ~ 10 9 K and ~ 5 x 10"K which confirms that cool, one temperature, pure Keplerian accretion solution is not unique. Indeed 
Eardley & Lightman (1975) found that a Keplerian disc is unstable due to thermal and viscous effects when viscosity parameter a (Shakura 
& Sunyaev 1973) is constant. Later Eggum et al. (1985) showed by numerical simulations that the Keplerian disc with a constant a collapses. 

Around eighties, therefore, the idea of two component accretion disc started floating around. For example, Paczynski & Wiita (1980) 
described a geometrically thick regime of the accretion disc in the optically thick limit, while Rees et al. (1982) introduced accretion torus in 
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the optically thin limit. Moreover the idea of sub-Keplerian, transonic accretion was introduced by Muchotrzeb & Paczyhski (1982), which 
was later improved by other authors (Chakrabarti 1989, 1996; Mukhopadhyay 2003). Other models were proposed by e.g. Gierlihski et al. 
(1999), Coppi (1999), Zdziarski et al. (2001), including a secondary component in the accretion disc. On the other hand, Narayan & Yi 
(1995) introduced a two temperature disc model in the regime of inefficient cooling resulting in a vertical thickening of the hot disc gas. 
Here the pressure forces are expected to become important in modifying the disc dynamics which is likely to be sub-Keplerian. Other models 
with similar properties were proposed by, e.g., Begelman (1978), Liang & Thompson (1980), Rees et al. (1982), Eggum, Coroniti & Katz 
(1988). Abramowicz et al. (1988) proposed a height-integrated disc model, namely "slim disc", having high optical depth of the accreting 
gas at super-Eddington accretion rate such that the diffusion time is longer than the viscous time. The model was further applied to study the 
thermal and viscous instabilities in optically thick accretion discs (Wallinder 1991; Chen & Taam 1993). 

Shapiro, Lightman & Eardley (1976) initiated a two temperature Keplerian accretion disc at a low mass accretion rate which is optically 
thin and significantly hotter than the single temperature Keplerian disc of Shakura & Sunyaev (1973). The optically thin hot gas cools down 
through the bremsstrahlung and inverse-Compton processes and could explain various states of Cyg X-l (Melia & Misra 1993). Similarly, 
the "ion torus" model by Rees et al. (1982) was applied to explain AGNs at a low mass accretion rate. However, the two temperature 
model solutions by Shapiro, Lightman & Eardley (1976) appear thermally unstable. Narayan & Popham (1993) and subsequently Narayan 
& Yi (1995) showed that introduction of advection may stabilize the system. However, the solutions of Narayan & Yi (1995), while of 
two temperatures, could explain only a particular class of hot systems with inefficient cooling mechanisms. They also described the hot 
flow based on the assumption of "self similarity" which is just a "plausible choice". They kept the electron heating decoupled from the 
disc hydrodynamical computations which merely is an assumption. Later on, the solutions were attempted to generalize by Nakamura et 
al. (1997), Manmoto et al. (1997), Medvedev & Narayan (2001), relaxing efficiency of cooling into the systems, but concentrating only 
on specific classes of solutions. On the other hand, Chakrabarti & Titarchuk (1995) and later Mandal & Chakrabarti (2005) modeled two 
temperature accretion flows around Schwarzschild black holes in the general "advective paradigm", emphasizing possible formation of shock 
and its consequences therein. However, they also did not include the effect of electron heating self-consistently into the hydrodynamical 
equation, and thus the hydrodynamical quantities do not get coupled to the rate of electron heating (see also Rajesh & Mukhopadhyay 2009). 

In the present paper, we model a selfconsistent accretion flow in the regime of two temperature transonic sub-Keplerian disc (see also 
Sinha, Rajesh & Mukhopadhyay 2009; a brief version of the present work, but around Schwarzschild black holes). We consider all the 
hydrodynamical equations of the disc along with thermal components and solve the coupled set of equations selfconsistently. We neither 
restrict to the advection dominated regime nor the self-similar solutions. We allow the disc to cool selfconsistently according to the thermo- 
hydrodynamical evolution and compute the corresponding cooling efficiency factor as a function of radial coordinate. We investigate that 
when does the disc switch from the radiatively inefficient nature to general advective paradigm and vice versa. 

In order to implement our model to explain observed sources, we focus on the under-luminous AGNs and quasars (e.g. Sgr A*), ultra- 
luminous quasars and highly luminous AGNs (e.g. PKS 0743-67) and ultra-luminous X-ray (ULX) sources (e.g. SS433), when the last 
items are likely to be the "radiation trapped" accretion discs. While the first two cases correspond to respectively sub-Eddington and super- 
Eddington accretion flows around supermassive black holes, the last case corresponds to super-Eddington accretors around stellar mass black 
holes. 

In the next section, we discuss the model equations describing the system and the procedure to solve them. Subsequently, we discuss the 
two temperature accretion disc flows around stellar mass and supermassive black holes, respectively in §3 and §4, for both sub-Eddington, 
Eddington and super-Eddington accretion rates. Section 5 compares the disc flow of low Shakura-Sunyaev (1973) a with that of high a and 
then between the flows around co and counter rotating black holes. Then we discuss the implications of the results with a summary in §6. 



2 MODEL EQUATIONS DESCRIBING THE SYSTEM AND SOLUTION PROCEDURE 

For the present purpose, we set five coupled differential equations describing the law of conservation in the sub-Keplerian optically thin 
accretion regime. Necessarily the set of equations describes the inner part of the accretion disc where the gravitational potential energy 
dominates over the centrifugal energy of the flow. 

Throughout, we express all the variables in dimensionless units, unless stated otherwise. The radial velocity i? and sound speed c s are 
expressed in units of light speed c, the specific angular momentum A in GM/c, where G is the Newton's gravitational constant and M is 
the mass of the compact object, for the present purpose black hole, expressed in units of solar mass M e , the radial coordinate x in units 
of GM/c 2 , the density p and the total pressure P accordingly. The disc fluid under consideration consists of ions and electrons — thus two 
fluid/temperature system, apart from radiation. Furthermore, at the high temperature, the disc flow with ions/electrons behaves as (almost) 
noninteracting gas. 



2.1 Conservation laws 

(a) Mass transfer: 
1 d 

-—(xp&) = 0, 
X ox 



(1) 
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whose integrated form gives the mass accretion rate 

M = -4nx!&, (2) 
where the surface density 

I = I„ph(x), (3) 

/„ = (Tn\fl(2n + 1)! (Matsumoto et al. 1984), (4) 

n is the polytropic index which is equal to l/(y- 1) when y is the ratio of specific heats, and half-thickness, based on the vertical equilibrium 
assumption, of the disc 

h(x) = c s x l/2 F~ 112 . (5) 

(b) Radial momentum balance: 
cW 1 dP A 2 

+ + F = (6) 

dx p ax x J 

when following pseudo-Newtonian approach of Mukhopadhyay (2002) 

p _ (x 2 -2a^x + a 2 ) 2 
x 3 [^(x-2) + a] 2 ' 

where a is the specific angular momentum (Kerr parameter) of the black hole. We also define a parameter 
gas pressure P ms 6y - 8 

j3 = — = — (e.g. Ghosh & Mukhopadhyay 2009), (8) 

total pressure P 3(y-l) 

where y may range from 4/3 to 5/3, P ga s = P t (ion pressure) + P e (electron pressure), such that 

*"£(£ + ^)='* *> 
/3c- \Hitni iieniij 

where T„ T e are respectively the ion and electron temperatures in Kelvin, m, is the mass of proton in gm, yu, and fi e respectively are the 
corresponding effective molecular weight, k the Boltzmann constant. We assume /? (and then y) constant throughout the flow. 

(c) Azimuthal momentum balance: 

"i - ss^-O- 

where following Mukhopadhyay & Ghosh (2003; hereinafter MG03) the shearing stress can be expressed in terms of the pressure and density 

as 

= -a(l„ +l P eq + I„& 2 p ell )h(x), (11) 

where a is the dimensionless viscosity parameter and P eq and p eq are the pressure and density respectively at the equatorial plane. Note that 
we will assume P eq ~ P and p eq ~ p in obtaining solutions. 

(d) Energy production rate: 

- 1 \ dx p dx J 
where following MG03 

Q + = a(/ n+1 P + I„& 2 p)h(x)^, (13) 
dx 

which is the heat generated by viscous dissipation, and Q le is the Coulomb coupling (Bisnovatyi-Kogan & Lovelace 2000) given in dimen- 
sionful unit as 

m^' 2 e\n e (T e T,\-"\ 

q le = — + — ln(A) (T, - T e ) erg/cm /sec^ 

m,m e \m e m,j 

when q ie = Q ie c n /(hG 4 M 3 ). (14) 

Here ;i, and n e denote number densities of ion and electron respectively, e the charge of an electron, ln(A) the Coulomb logarithm. We also 
define (MG03) 

T l - B 

r3 = l + 4-Tp> (15) 
= (4-3^(r-i) 

H P+ 12(7 - D(l -0) 
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(e) Energy radiation rate: 

m*{dp. Tl ^) = Qk . Q -, (17) 



— \ \ dx p dxj 

where QT is the heat radiated away by the bremsstrahlung {qt r ), synchrotron (q sy „) processes and inverse Comptonization (q comp ) due to soft 
synchrotron photons, given in dimensionful form as 

q~ = q br + + q comp , when q' = Q~ c ll /(h G 4 M 3 ). (18) 

Various components of the cooling processes may be read as (see Narayan & Yi 1995; Mandal & Chakrabarti 2005 for detailed description, 
what we do not repeat here) 

qtr = 1.4 x 10~ 27 n e ?nT ,/2 (1 + 4.4 x 10" 10 r f ) erg/cm 3 /sec, 
In v 3 

9.1.™ = 3^2^" "/j er g/ cm3 / sec ' ^ = xGM/c 2 , 

qcomp = Tq s ,„, T = 1]i II - I — I I, m = (1 _ A y P = 1 - exp(-T es ), 

ln(A) 

where T es is the scattering optical depth given by 

T es = K es p h (20) 

where k„ = 0.38 cm 2 /gm and v a is the synchrotron self-absorption cut off frequency determined by following Narayan & Yi (1995). Note that 
without a satisfactory knowledge of the magnetic field in accretion disks, following Mandal & Chakrabarti (1995), we assume the maximum 
possible magnetic energy density to be the gravitational energy density of the flow. As the total optical depth should include the effects of 
absorption due to nonthermal processes, effective optical depth is computed as 



A = 1 + 40 e + 166*;, e = kT e /m e c 2 , 772 = — ( 1 + , x a = hvjm e c 2 , (19) 



r eff — V r fJ T abs (21) 

where 

h , \ GM 

T abs ~ ^ a j4 {4t>r + Isyn + qcomp J (22) 

when <r is Stefan-Boltzmann constant. 

Now combining all the above equations we obtain 

d& N(x,&,c s ,A,T e ) 



dx D(&, c s ) 

where 

N = -&~c s J H\—c s + & - a — — 2HJ + tvct s L + aH\- 



(23) 



r 3 - i " x \ /„ 1 / /„ r 3 - l 1 \ x 



+ ^#»" 2 , (24) 
M 



and 

L = (s - ■ H = M + ■ 7 = ( C " L + 7 - F ) ■ (26) 

We know that around the sonic radius N = D = (Mukhopadhyay 2003) and hence obtain the Mach number at this critical radius from 
D = 



where 

a T[ + 1 , 2(n + 1) 2 /„+i 

A = 2a (7„ +1 - /„) , B = — 2a — (7„ + i - /„) , 

1 - 1 3 1 - 1 3 /„ 

C = a 2 ^i(l - 2/„ +l ). (28) 
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Also from N = 0, we can compute explicitly c sc as a function of sonic/critical radius x c . For a physical x c , what one has to adjust in order 
to obtain a physical solution connecting outer boundary to black hole horizon through x c , c sc and then & c can be assigned, discussed in 
APPENDIX A in detail. Note that an improper x c may lead to an unphysical/imaginary c sc and & c . 
Finally combining Eqns. {6^, J 1 Ob and | |17I > we obtain 

dc s _ lc s &\d& J 
dx \& c, j dx c s ' 

dX (2axl„ +l (c] \ \d& !c]-2xaJ \ 

■&c s \ + ax — + — +■& , (30) 



dx \ i?c s 1„ \& ] ] dx \ c 

f = d- + (i - W (4 + l) + - Q-) ■ (3D 

dx cj dx \c- J M F 1 ' 2 v ' 

Hence knowing & we can obtain other variables c s , A, T e . Note that d&/dx is indeterminate (of 0/0 form) at x c . APPENDIX B discusses the 
procedure to obtain d&/dx at x c . 

As the entropy increases inwards in advective flows (e.g. Narayan & Yi 1994; Chakrabarti 1996; MG03), there is a possibility of 
convective instability and then corresponding transport, as proposed by Narayan & Yi (1994). Dynamical convective instability arises when 
the square of effective frequency 

^ ff = v|v + < (32) 
such that y BV is the Brunt-Vaisala frequency given by 

BV p dx dx ( p ) 
and v r is the radial epicyclic frequency. 



2.2 Solution procedure 

In order to obtain the steady state solution, as of previous work (MG03, Mukhopadhyay 2003), primarily we need to fix the appropriate 
critical radius x c (in fact the energy at the critical radius which is not conserved in the present cases) and the corresponding specific angular 
momentum A c of the flow. The detailed description of the procedure to obtain physically meaningful values of x c and A c , to be determined 
iteratively, is given in APPENDIX A. As the flow is considered to be of two temperatures, at x c an appropriate electron temperature T ec also 
needs to be determined; also discussed in APPENDIX A. Note that one has to adjust the set of values x c , A c , T ec appropriately/iteratively to 
obtain self-consistent solution connecting outer boundary and black hole event horizon through x c . Depending on the type of accreting system 
to model, we then have to specify the related inputs: M, M, y and a. Important point to note is that unlike former works (e.g. Chakrabarti 
& Titarchuk 1995, Chakrabarti 1996, MG03) here x c changes with the change of M, because the various cooling processes considered here 
explicitly depend on M. Finally, we have to solve the Eqn. i23l from x c to inwards — upto the black hole event horizon, and to outwards 
— upto the transition radius x where the disc deviates from the Keplerian to the sub-Keplerian regime such that A/A K = 1 (A K being the 
specific angular momentum of the Keplerian part of the disc). FigureQ]shows h° w the ratio A/A K varies as a function of radial coordinate for 
different a. Note that higher a, which corresponds to a lower disc angular momentum (Mukhopadhyay 2003), reassembles the Keplerian part 
advancing with a smaller size of the sub-Keplerian disc. On the other hand, for a lower a the inner edge of the Keplerian component recedes. 
The fact of moving in and out of the inner edge of the disc reassembles respectively the soft and hard state of the black hole (e.g. Gilfanov et 
al. 1997). Hence it is naturally expected to link with the spin of the black hole. 

However, the important point to note is that there is no selfconsistent model to describe the transition region where A/A K = 1. Therefore, 
the transition of the flow from the Keplerian to sub-Keplerian regime does not appear smooth. This is mainly because the set of equations 
used to model the sub-Keplerian flow is not valid to explain the cold Keplerian flow, unless an extra boundary condition is imposed at the 
outer edge of the sub-Keplerian disc. However, in the present paper we do not intend to address the transition zone; rather we prefer to stick 
with the sub-Keplerian flow. Narayan et al, (1997) imposed boundary conditions at both the ends of accretion flows along with at the critical 
radius to fix the problem, at the cost of more input parameters than the parameters chosen in the present work. But still the transition of the 
flow from the Keplerian to sub-Keplerian regime remains undefined. Later Yuan (1999) discussed how the solutions vary with the change of 
outer boundary conditions influencing the structure of an optically thin accretion flow. 

Below we discuss solutions in various parameter regimes to understand properties of the accretion disc around, first, stellar mass 
(M = 10) and then super-massive (M = 10 7 ) black holes. 
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log(x) 

Figure 1. Variation of ratio of disc specific angular momentum to corresponding Keplerian angular momentum as a function of radial coordinate, when solid, 
dashed, dotted, dot-dashed lines correspond to the cases with a = -0.5, 0, 0.5, 0.998 respectively. Other parameters are M = 1, M — 10, a = 0.01. 

3 TWO TEMPERATURE ACCRETION DISC AROUND STELLAR MASS BLACK HOLES 

Primarily we concentrate on two extreme regimes: (1) sub-Eddington and Eddington limits of accretion, (2) super-Eddington accretion. 
Furthermore, at each case of accretion rate, we focus on solutions around nonrotating (Schwarzschild) and rotating (Kerr with a = 0.998) 
black holes. 

One of our aims is to understand how the explicit cooling processes affect the disc dynamics and then the cooling efficiency vary over 
the disc radii. The cooling efficiency / is defined as the ratio of the energy advected by the flow to the energy dissipated, which is 1 for the 
advection dominated accretion flow (in short ADAF; Narayan & Yi 1994, 1995) and less than 1 for the general advective accretion flow (in 
short GAAF; Chakrabarti 1996; Mukhopadhyay 2003; MG03) in general. Therefore, / directly controls the ion and electron temperatures of 
the disc. Far away from the black hole where the gravitational power is weaker, the angular momentum profile becomes Keplerian and thus 
the disc becomes (or tends to become) of one temperature in the presence of efficient cooling. 

3.1 Sub-Eddington and Eddington accretors 

3.1.1 Schwarzschild black holes 

We first consider flows around static black holes where the Kerr parameter a = 0. Figure [2] shows the behavior of flow variables as functions 
of radial coordinate for M = 0.01,0.1, 1; throughout in the text we express M in units of Eddington limit. The sets of input parameters for 
the model cases described here are given in Table 1. Figured 
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Figure 2. Variation of dimensionless (a) radial velocity, (b) density, (c) cooling factor, (d) square of convective frequency, as functions of radial coordinate 
for sub-Eddington and Eddington accretion flows. Solid, dashed, dotted curves are for M = 0.01,0.1, 1 respectively. Other parameters are a = 0, a = 0.01, 
M = 10; see Table 1 for details. 



Table 1: Parameters for accretion with a = 0.01 around black holes of M = 10, when the subscript 'c' indicates the quantity at the critical 
radius and T ec is expressed in units of m,c 2 /£ 



M 


a 


r 


X c 




Tec 






Sub-Eddington, 


Eddington 


accretors 




0.01 





1.5 


5.5 


3.2 


0.0001 


0.01 


0.998 


1.5 


3.5 


1.7 


0.0001 


0.1 





1.4 


5.5 


3.2 


0.000164 


0.1 


0.998 


1.4 


3.5 


1.7 


0.000153 


1 





1.35 


5.5 


3.2 


0.000225 


1 


0.998 


1.35 


3.5 


1.7 


0.0002122 






Super-Eddington 


accretors 






10 





1.345 


5.5 


3.2 


0.000181565 


10 


0.998 


1.345 


3.5 


1.7 


0.0004432 


100 





1.34 


5.5 


3.2 


0.00038678 


100 


0.998 


1.34 


3.5 


1.7 


0.00055 



verifies that a higher radial velocity corresponds to a lower mass accretion rate of the flow (~ 0.01) which results in less possible 
accumulation of matter in a particular radius attributing to a lower disc density (Fig. EJs). This hinders the bremsstralung process to cool the 
flow. However, at around x = 30 the centrifugal barrier dominates and brings the velocity & down, particularly for M = 0.01, which finally 
merges with that of higher Ms. On the other hand, a lower M corresponds to a gas dominated hot flow, which is radiatively less efficient and 
quasi-spherical in nature. As a result j? is high, as seen in Fig. [2^. Efficiency of cooling is shown in Fig.[2j;. Naturally a low M corresponds 
to a radiatively inefficient flow rendering f>0.9 upto x = 30. At x < 30, the dominance of centrifugal barrier slows down the infall which 
increases the residence time of matter in the disc before plunging into the black hole. This allows matter to have enough time to radiate by the 
synchrotron process and inverse Comptonization due to synchrotron soft photons, rendering / — > close to the black hole. In other words, 
for M = 0.01, the disc is essentially radiatively inefficient, upto x ~ 30, and therefore the electron temperature never goes down. However, 
the density sharply increases in the vicinity of the black hole (Fig.[2p) which favors efficient cooling at a high temperature. 

Therefore, although far away from the black hole a sub-Eddington flow appears to be radiatively inefficient, from x = 30 onwards it 
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(a) (b) 




0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 

log(x) log(x) 



Figure 3. Variation of (a) dimensionless energy of Coulomb coupling (thicker line), bremsstrahlung (dotted line), synchrotron (solid), inverse Comptonization 
due to synchrotron photon (dashed line) processes in logarithmic scale, (b) corresponding ion (solid) and electron (dotted) temperatures in units of 10 8 K, as 
functions of radial coordinate for M = 0.01. (c), (e) Same as (a) except M = 0.1, 1 respectively, (d), (f) Same as (b) except M = 0.1, 1 respectively. Other 
parameters are a = 0, a = 0.01, M = 10; see Table 1 for details. 

turns out to be a radiatively efficient advective flow with / much less than unity. However, for M = 0. 1 , 1 the density is higher than that for 
M = 0.01 and hence the bremsstrahlung effect starts playing role in radiation mechanisms at much outer radii. This decreases the ion-electron 
temperature difference at the transition radius x . However, as the flow advances the synchrotron and corresponding inverse Compton effects 
dominate attributing to strong radiative loss. This renders /J0.5 upto x = 10. Further in, a strong radial infall, in absence of any centrifugal 
barrier compared to a low M case, does not permit matter to radiate enough, rendering / upto 1. This is particularly because the strong 
advection decreases the residence time of the flow before plunging into the black hole and thus renders a weaker ion-electron coupling. This 
in turn hinders the transfer of energy from the ions to electrons attributing ions to remain hot, while electrons continue to be cooled down 
further by radiative processes. 

However, Fig.|2jl shows that either of the cases do not exhibit convective instability (see, however, Narayan, Igumenshchev & Abramow- 
icz 2000, Quataert & Gruzinov 2000) upto very inner edge, evenif the radiatively inefficient flow deviates to a radiatively efficient GAAF. At 
a very inner edge, discs with M = 0.1, 1 particularly appear to be marginally unstable, which, although, seems not playing any role in angular 
momentum transfer. 

Figure[3]describes how do the various cooling processes and corresponding temperature profiles vary as functions of radial coordinate. 
At a low M (= 0.01) the system is radiatively inefficient, relative to that of higher Ms (= 0.1, 1), which brings out a hot two temperature 
Keplerian-sub-Keplerian transition region. As the flow advances through the sub-Keplerian part, strong two temperature nature remains 
intact. For a higher M (= 1), however, the transition region is of marginally two temperatures. This is due to the efficient bremsstrahlung 
radiation at high density. In the vicinity of black hole T e ~ 10 9 K in the flow with M = 0.01 when / — > 0, as explained above, in the contrary 
to the cases with M = 0.1, 1 when / — > 1 and T e sharply decreases. Note that the accretion disc around a stellar mass black hole is arrested 
by significant magnetic field. This results in dominance of the synchrotron effect over the bremsstrahlung as the flow advances. 

3.1.2 Kerr black holes 

We consider the rotating black holes with a = 0.998. As discussed earlier (Mukhopadhyay 2003), the angular momentum of the flow should 
be smaller around a rotating black hole compared to that around a static black hole. This reassembles advancing the Keplerian component 
which decreases size of the sub-Keplerian part. As the disc remains Keplerian (which is radiatively efficient) upto, e.g., x ~ 100 (see the 
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outer radius in Fig. [4j», the flow cools down significantly before deviating to the sub-Keplerian regime. Figure [4^ shows that the velocity 
profiles for all M-s are similar to each other, in absence of a strong centrifugal barrier. However, /, while very small at x ~ 100, increases as 
the sub-Keplerian flow advances. This is because the residence time of the flow decreases in an element of the sub-Keplerian disc hindering 
cooling processes to complete. Figure^ shows that even / — » 0.9 at x ~ 30 for M = 0.01, when the density is lowest (see Fig.|4p). 

However, as the flow approaches to the black hole the synchrotron emission increases, and hence the system acquires enough soft 
photons which help in occurring the inverse Compton process. As a result the flow cools down further. When M = 0.01 the cooling process 
at the very inner edge of the accretion disc is dominant due to relatively high residence time of the flow, compared to that of a higher M, 
rendering high T e and then / — > very close to the black hole. Figure |4jl proves that the flow is convectionally stable all the way upto the 
black hole event horizon. For MJ0.1, on the other hand, the flow is strongly advective and then unable to cool down before plunging into the 
black hole. Figure [5] shows the profiles of cooling processes and ion-electron temperatures. Basic nature of the profiles is pretty similar to 
that of Schwarzschild cases, except all of them advance in. 



3.2 Super-Eddington accretors 

The "radiation trapped" accretion disc can be attributed to the radiatively driven outflow or jet. This is likely to occur when the accretion 
rate is super-Eddington (Lovelace et al. 1994, Begelman et al. 2006, Fabbiano 2004, Ghosh & Mukhopadhyay 2009), as seen in the ultra 
luminous X-ray (ULX) sources such as SS433 (with luminosity ~ 10 40 erg/s or so; Fabrika 2004). In order to describe such sources, the 
models described below are the meaningful candidates. We consider M = 10, 100. 



3.2. 1 Schwarzschild black holes 

A high mass accretion rate significantly enhances density, upto two orders of magnitude compared to that of a low M, which severely affects 
/ and finally temperature profiles. The profiles of velocity shown in Fig.|6]are quite similar to that of sub-Eddington and Eddington cases. 
Because of similar reasons explained in §3.1.1 the profile exhibits a stronger centrifugal barrier for a lower M (= 10). A lower M flow 
will have relatively more gas and then quasi-spherical structure compared to that of a higher M (= 100), which results in a lower velocity 
in the latter case. On the other hand, a lower velocity corresponds to a higher density which results in the strong bremsstrahlung radiation 
rendering / — > 0. For a lower M, at x ~ 50, the energy radiated due to bremsstrahlung process becomes weaker than the energy transferred 
from protons to electrons through the Coulomb coupling (see Fig. [7}, which increases / (see Fig. [6};). Subsequently, the synchrotron process 
becomes dominant (see Fig. [7), reassembling / — > 0. However, very close to the black hole a strong advection does not allow the flow, 
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independent of M, to radiate efficiently rendering / — » 1 again. This also results in marginal convective instability at x < 10, as shown in 
Fig. |6jl. 

Figure[7]shows that discs remain of one temperature at the transition radius. As the flows advance with a sub-Keplerian angular momen- 
tum, T p profile deviates from that of T e . For M = 100, the high density flow is dominated by very efficient bremsstrahlung radiation all the 
way. In the vicinity of the black hole, an efficient cooling reassembles a sharp downfall of T e . 



3.2.2 Kerr black holes 

As A decreases in the case of a higher a, like low M cases (see Mukhopadhyay 2003), the transition region advances an order of magnitude 
compared to that of Schwarzschild black holes. Similar to the cases of low M, as shown in Fig. [8ji, any centrifugal barrier smears out. 
However, unlike the flow around a static black hole, here the disc with M = 10 remains stable upto very close to the black hole. The reason is 
that a high a corresponds to a larger inner edge of the disc and thus the residence time of matter in the disc is higher. As a result the radiative 
processes keep cooling and then stabilizing the flow upto very inner edge. 

Figure [9] shows that although a high M exhibits a one temperature transition zone due to extremely efficient cooling processes, par- 
ticularly due to bremsstrahlung radiation, as M decreases the Keplerian disc itself becomes of two temperatures before deviating to the 
sub-Keplerian zone, unlike that of the Schwarzschild case. This is mainly because a flow with a high a brings the Keplerian disc further in 
where the transport of angular momentum increases leading to the decrease of the residence time of matter which does not allow an efficient 
cooling. However, the basic behaviours of various cooling processes is pretty similar to that around a static black hole. 



4 TWO TEMPERATURE ACCRETION DISC AROUND SUPERMASSIVE BLACK HOLES 

As of stellar mass black holes, here also we concentrate on two extreme regimes: (1) sub-Eddington and Eddington limits of accretion, (2) 
super-Eddington accretion; focusing on both nonrotating (Schwarzschild) and rotating (Kerr with a = 0.998) black holes. 
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(a) (b) 




log(x) log(x) 

Figure 6. Variation of dimensionless (a) radial velocity, (b) density, (c) cooling factor, (d) square of convective frequency, as functions of radial coordinate for 
super-Eddington accretion flows. Solid, dashed curves are for M = 10, 100 respectively. Other parameters are a = 0, a = 0.01, M = 10; see Table 1 for details. 



(a) (b) 




1 1.5 2 1 1.5 2 

log(x) log(x) 

Figure 7. Variation of (a) dimensionless energy of Coulomb coupling (thicker line), bremsstrahlung (dotted line), synchrotron (solid), inverse Comptonization 
due to synchrotron photon (dashed line) processes in logarithmic scale, (b) corresponding ion (solid) and electron (dotted) temperatures in units of 10 8 K, as 
functions of radial coordinate for M = 10. (c) Same as (a) except M = 100. (d) Same as (b) except M = 100. Other parameters are a = 0, a = 0.01, M = 10; 
see Table 1 for details. 




Figure 9. Same as Fig.|7] except a = 0.998. 
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log(x) log(x) 

Figure 10. Variation of dimensionless (a) radial velocity, (b) density, (c) cooling factor, (d) square of convective frequency, as functions of radial coordinate 
for sub-Eddington and Eddington accretion flows. Solid, dashed, dotted curves are for M = 0.01,0.1, 1 respectively. Other parameters are a = 0, a = 0.01, 
M = 10 7 ; see Table 2 for details. 

4.1 Sub-Eddington and Eddington accretors 

The under-luminous AGNs and quasars (e.g. Sgr A*) had been already described by advection dominated model, where the flow is expected 
to be substantially sub-critical/sub-Eddington with a very low luminosity (JlO 35 erg's). Therefore the present cases, particularly of M;0.01, 
could be potential models in order to describe under-luminous sources. 

4.1.1 Schwarzschild black holes 

Table 2 lists the sets of input parameters for the model cases described here. Naturally a disc around a supermassive black hole will have 
much lower density compared to that around a stellar mass black hole. Therefore, the cooling processes, particularly the bremsstrahlung 
radiation which is density dependent, are expected to be inefficient leading to a high /. However, the velocity profiles shown in Fig.llQfa are 
very similar/same to that around a stellar mass black hole. Figure [TOb shows that / — > 1 in most of the sub-Keplerian regime for M = 0.01. 
As M increases, the density increases and thus the bremsstrahlung radiation increases, as shown in Fig.[TT] which leads to the transition 
of radiatively inefficient flow to GAAF. When M = 1 the bremsstrahlung effect is very high resulting an GAAF with / much smaller than 
unity upto very close to the black hole. Figure [Toti shows that close to the black hole there is a possible convective instability for all M-s. 
This is because a strong advection of matter close to the black hole hindering cooling processes which results in / — > 1. This reassembles a 
possible convective instability at the inner edge. For a higher M, the density is high which favours convection and thus brings the convective 
instability earlier in, at a relatively outer radius. 
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Figure 11. Variation of (a) dimensionless energy of Coulomb coupling (thicker line), bremsstrahlung (dotted line), synchrotron (solid), inverse Comptonization 
due to synchrotron photon (dashed line) processes in logarithmic scale, (b) corresponding ion (solid) and electron (dotted) temperatures in units of 10 8 K, as 
functions of radial coordinate for M = 0.01. (c), (e) Same as (a) except M = 0.1, 1 respectively, (d), (f) Same as (b) except M = 0.1, 1 respectively. Other 
parameters are a = 0, a = 0.01, M = 10 7 ; see Table 2 for details. 



Table 2: Parameters for accretion with a = 0.01 around black holes of M = 10 7 , when the subscript 'c' indicates the quantity at the critical 
radius and T ec is expressed in units of m,c 2 /£ 



M 


a 


y 


x c 


A c 








Sub-Eddington, 


Eddington 


accretors 




0.01 





1.5 


5.5 


3.2 


0.0001 


0.01 


0.998 


1.5 


3.5 


1.7 


0.0001 


0.1 





1.4 


5.5 


3.2 


0.000178 


0.1 


0.998 


1.4 


3.5 


1.7 


0.00023 


1 





1.35 


5.5 


3.2 


0.0002493 


1 


0.998 


1.35 


3.5 


1.7 


0.000295 






Super-Eddington 


accretors 






10 





1.345 


5.5 


3.2 


0.000427 


10 


0.998 


1.345 


3.5 


1.7 


0.0006 


100 





1.34 


5.5 


3.2 


0.0003874 


100 


0.998 


1.34 


3.5 


1.7 


0.00059 



The temperature profiles shown in Fig. [TT] are pretty similar to what we obtain in stellar mass black holes what we do not explain 
here again. However, note that unlike stellar mass black holes, only the bremsstrahlung radiation is effective in cooling the flow around a 
supermassive black hole, particularly for M = 1 . 



4. 1.2 Kerr black holes 

The specific angular momentum of the black hole is chosen to be a = 0.998. The basic hydrodynamical properties are similar to that in 
Schwarzschild cases, except, like the flows around stellar mass black holes, the transition region advances due to a smaller angular momentum 
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Figure 12. Same as Fig.[TS] except a = 0.998. 



of the flow, shown in Fig. [12] As discussed in §3.2.2, a higher a corresponds to a smaller A which in turn decreases j? at a particular radius 
of the inner edge of the disc, when the inner edge is stretched in, compared to that around a Schwarzschild black hole. This results in the 
increase of residence time of the flow in the sub-Keplerian disc before plunging into the black hole. Therefore, the bremsstrahlung process 
keeps cooling and then stabilizing the flow upto very inner edge. Important point to note, as a consequence, is that the disc flow around a 
rotating black hole is convectively more stable compared to that around a nonrotating black hole. This is particularly because the density 
gradient of the inner (e.g. the vicinity of x = 2, which is the event horizon for a nonrotating black hole) flow around a rotating black hole is 
less stepper compared to that around a nonrotating black hole, and hence the flow is convectively more stable in the former case. 

Figure [J_3] shows that basic features of the temperature profiles are similar to the cases of static black holes. However, the transition 
region reveals that for a lower M, the Keplerian flow exhibits the inverse Comptonization via synchrotron photons. As the flow advances with 
a sub-Keplerian angular momentum, the residence time decreases and thus inverse Comptonization decreases, resulting a hotter flow. 



4.2 Super-Eddington accretors 

Ultra-luminous accretors with a high kinetic luminosity (~ 10 46 - 10 49 erg/s) radio jet have been observed in the highly luminous AGNs 
and ultra-luminous quasars (e.g. PKS 0743-67; Punsly & Tingay 2005), possibly in ULIRs (Genzel et al. 1998) and narrow-line Seyfert 1 
galaxies (e.g. Mineshige et al. 2000). Therefore, the following cases could be potential models to explain such sources. 



4.2. 1 Schwarzschild black holes 

Figures [74l and [731 show that the basic flow properties are pretty similar to that around stellar mass black holes, except in the present cases 
the centrifugal barrier smears out. This is because a high black hole mass corresponds to a low density of the flow and thus a fast infall. This 
also results, unlike stellar mass black holes, in an inefficient synchrotron radiation even at the inner edge of the disc. 



4.2.2 Kerr black holes 

Figures[J_6]and[T7]repeat the same story of that around stellar mass black holes, but with the smeared centrifugal barrier, as described above 
for static black holes. However, due to decreasing density, the overall cooling effects decrease keeping the disc hotter, particularly for M = 10. 
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5 COMPARISON BETWEEN FLOWS WITH DIFFERENT (X AND AROUND COROTATING AND COUNTER ROTATING 
BLACK HOLES 

So far we have restricted to a typical Shakura-Sunyaev viscosity parameter a = 0.01, for corotating black holes. Now we plan to explore a 
lower a , as well as a counter rotating black hole to understand any significant change in the flow behaviour. 



5.1 Comparison between flows with a = 0.01 and a = 0.0001 

Decreasing a naturally decreases the rate of energy-momentum transfer between any two successive layers of the fluid element and increasing 
the residence time of the flow in the sub-Keplerian disc. This also recedes the Keplerian-sub-Keplerian transition region further out. This is 
mainly because a low value of a can not keep continuing the outward angular momentum transport efficiently in the Keplerian flow below 
a certain radial coordinate. Therefore, the disc flow can not remain Keplerian and becomes sub-Keplerian at a larger radius, compared to a 
flow of high a. 

We know, on the other hand, that increasing residence time increases the possibility of completing various radiative processes in the disc 
flow, before the infalling matter plunges into the black hole. Therefore, the flow is expected to appear cooler with smaller /. Hence, for the 
purpose of comparison, we consider a flow with M = 0.01 around a supermassive black hole of M = 10 7 , e.g. Sgr A", which is radiatively 
inefficient and hot for a = 0.01. 

Figure Q~8] shows that although the velocity profiles are similar for both the values of a, the size of the sub-Keplerian disc is about five 
times for a = 0.0001 than that for a = 0.01. Inside x = 17 the low a disc flow becomes cooler very fast, rendering / — > at x > 10 (see Fig. 
I18fc). Therefore, the flow sharply transits from radiatively inefficient in nature to GAAF. As a consequence, the low a flow remains stable, 
as shown in Fig. 1 18b. all the way upto the event horizon. As the sub-Keplerian flow of a smaller a extends further away where the influence 
of black hole is very weak, T e and T, merge (see Fig.ll8U) before the flow crosses the transition radius, unlike the flow with a = 0.01 when 
Tj > T e there. 
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Figure 14. Variation of dimensionless (a) radial velocity, (b) density, (c) cooling factor, (d) square of convective frequency, as functions of radial coordinate 
for super-Eddington accretion flows. Solid, dashed curves are for M = 10, 100 respectively. Other parameters are a = 0, a = 0.01, M = 10 7 ; see Table 2 for 
details. 



5.2 Comparison between flows around co and counter rotating black holes 

As we already discussed that the model with a low mass accretion rate around a supermassive black hole is a potential case in explaining the 
observed dim source Sgr A*. On the other hand, the ultra-luminous X-ray sources presumably correspond to the models with a high mass 
accretion flow around a stellar mass black hole. Therefore, in order to compare the flow properties between co and counter rotating black 
holes, we choose these two extreme cases. 

Qualitatively, the flows with similar initial conditions around co and counter rotating black holes of same mass are similar, as shown 
in Figs. [19] EfJ] for a = ±0.5. However, the sub-Keplerian disc size around the black hole with a = -0.5 is smaller due to smaller value 
of the effective angular momentum (Mukhopadhyay 2003) of the system. Hence the radial velocity is almost an order of magnitude higher, 
particularly at the inner edge, for a = -0.5. 



6 DISCUSSION AND SUMMARY 

We model the two temperature accretion flow, particularly around black holes, combining the equations of conservation and comprehensive 
cooling processes. We consider self-consistently the important cooling mechanisms: bremsstrahlung, synchrotron and inverse Comptoniza- 
tion due to synchrotron photons, where ions and electrons are allowed to have different temperatures. As matter falls in, hot electrons cool 
through the various cooling mechanisms, particularly by the synchrotron emission when the magnetic field is high. This is particularly the 
case for the flow around stellar mass black holes where the magnetic field may also act as a boost in transporting the angular momentum. 
However, in the present paper, we do not consider such processes in detail, rather stick with the standard a-prescription. 

By solving a complete set of disc equations we show that in general the disc system exhibits GAAF. However, in certain circumstances 
GAAF becomes radiatively inefficient, depending on the flow parameters and hence efficiency of cooling mechanisms. Transitions from 
GAAF to radiatively inefficient flow and vice versa are clearly explained by the cooling efficiency factor /, shown in each model cases. 
While the previous authors, who proposed ADAF (Narayan & Yi 1994, 1995), especially restricted with flows having / = 1 (inefficient 
cooling), here we do not impose any restriction to the flow parameters to start with and let the parameter / to determine self-consistently 
as the system evolves. Therefore, our model is very general whose special case may be understood as a radiatively inefficient advection 
dominated flow. 

We have explored especially the optically thin flows incorporating bremsstrahlung, synchrotron and inverse Comptonization processes. 
In Fig.[2T]we show the variation of the effective optical depth as a function of disc radii for two limiting cases. While flows around rotating 
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Figure 15. Variation of (a) dimensionless energy of Coulomb coupling (thicker line), bremsstrahlung (dotted line), synchrotron (solid), inverse Comptonization 
due to synchrotron photon (dashed line) processes in logarithmic scale, (b) corresponding ion (solid) and electron (dotted) temperatures in units of 10 8 K, as 
functions of radial coordinate for M = 10. (c) Same as (a) except M = 100. (d) Same as (b) except M = 100. Other parameters are a = 0, a = 0.01, M = 10 7 ; 
see Table 2 for details. 
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Figure 16. Same as Fig.[l4] except a = 0.998. 
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Figure 18. Comparison between solutions for high and low a: Variation of (a) velocity, (b) square of convective frequency, (c) cooling factor, (d) ion (upper set 
of lines) and electron (lower set of lines) temperatures, a functions of radial coordinate, when solid lines correspond to a = 0.01 and dashed lines correspond 
to a = 0.0001. Other parameters are M = 0.01, M = 10 7 , a = 0. 
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Figure 19. Comparison between solutions for co and counter rotating stellar mass black holes: Variation of (a) velocity, (b) square of convective frequency, 
(c) cooling factor, (d) ion (upper set of lines) and electron (lower set of lines) temperatures, as functions of radial coordinate, when solid lines correspond to 
b = 0.5 and dashed lines correspond to b = -0.5. Other parameters are fit = 100, M = 10, a = 0.01. 
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Figure 20. Same as Fig. [TO] except M = 0.01, M = 10 7 . 
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Figure 21. Variation of the effective optical depth as a function of radial coordinate. Solid (a = 0) and dotted (a = 0.998) curves correspond to M = 10, 
M = 100 and dashed (a = 0) and dot-dashed (a = 0.998) curves correspond to M = 10 7 , M = 0.01. See Tables 1 and 2 for details. 

black holes appear slightly thinner compared to the corresponding cases of static black holes, in general T c ffj5 x 1CT 4 . This verifies our choice 
of optically thin flows throughout. However, for the present purpose, when the main aim is to understand disk dynamics in the global, viscous, 
two-temperature regime, we have ignored inverse Comptonization due to bremssstrahlung photons, if any. This may be important in cases of 
very super-Eddington accretion flows which we plan to explore in future, particularly, in analysing the underlying spectra. 

The temperature of the flow depends on the accretion rate. If the accretion rate is low and thus the flow is radiatively inefficient, then the 
disc is hot. Such a hot flow is being attempted to model since 1976 (Shapiro, Lightman & Eardley 1976) when it was assumed that locally 
Q + ~ <2 and thus / — > 0. While the model was successful in explaining observed hard X-rays from Cyg X-l, it turned out to be thermally 
unstable. Rees et al. (1982) proposed a hot ion torus model avoiding / to unity. In the similar spirit Narayan & Yi (1995) proposed the hot 
two temperature solution in the assumption of / — > 1 including the strong advection into the flow. Abramowicz et al. (1995), based on the 
single temperature model, showed that the optically thin disc flow of accretion rate more than one Eddington does not have an equilibrium 
solution. However, they did not attempt to solve the complete set of differential equations. Based on some simplistic assumptions they showed 
the importance of advective cooling. Moreover, a single temperature description does not allow them to include all the underlying physics 
necessary to describe the cooling processes. In the due course, Mandal & Chakrabarti (2005) proposed a two temperature disc solution where 
the ion temperature could be as high as ~ 10 12 K. However, they particularly emphasized on how does the shock in the disc flow enable 
cooling through the synchrotron mechanism, without carrying out a complete analysis of the dynamics. The present paper describes, to 
our knowledge, the first comprehensive work to model the two temperature accretion flow self-consistently by solving the complete set of 
underlying equations without any pre-assumptive choice of the flow variables to start with. 

The generality lies not only in its construction but also its ability to explain the under-luminous to ultra-luminous sources, stellar mass 
to supermassive black holes. Table 3 lists the luminosities for a wide range of parameter sets, obtained by our model. It reveals that for a very 
low mass accretion rate M = 0.0001 around a supermassive black hole, the luminosity comes out to be L ~ 10 34 erg/sec, which indeed is 
similar to the observed luminosity from a under-luminous source Sgr A*. In other extreme, for M = 100 around a similar black hole, L ~ 10 47 
erg/sec, similar to what observed from the highly luminous AGNs like PKS 0743-67. On the other hand, when the black hole is considered to 
be of stellar mass, then at a high M = 100, the model reveals L ~ 10 40 erg/sec which is similar to the observed luminosity from ULX sources 
(e.g. SS433). 
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Table 3: Luminosity in erg/sec 



M 


M 


y 


L 


0.0001 


10 7 


1.6 


10 34 


0.01 


10 7 


1.5 


10 38 


1 


10 7 


1.35 


5 x 10 42 


100 


10 7 


1.34 


10 47 


0.01 


10 


1.5 


10 33 


1 


10 


1.35 


7 x 10 36 


100 


10 


1.34 


10 40 



In general, an increase of accretion rate increases density of the flow which may lead to a high rate of cooling and thus decrease of 
the cooling factor /. Hence, / is higher, close to unity reassembling radiatively inefficient flows, for sub-Eddington accretors, and is lower, 
sometimes close to zero, for super-Eddington flows. Actual value of / in a flow also depends on the behaviour of hydrodynamic variables 
which determine the rate of cooling processes. Naturally, as the flow advances from the transition region to the event horizon, / varies 
between and 1. However, if the black hole is considered to be rotating, the flow angular momentum decreases and thus the radial velocity 
increases. This in turn reduces the residence time of the sub-Keplerian flow hindering cooling processes to complete. This flow is then 
expected to be hotter and hence / to be higher compared to that around a static black hole. Therefore, the system may tend to be radiatively 
inefficient, evenif its counter part around a static black hole appears to be an GAAF. However, this also depends on the value of a, as shown 
in Fig. [THJ A low value of a increases the residence time of matter in the disc which helps in cooling processes to complete, rendering a 
radiatively inefficient flow to switch over to GAAF. This feature may help in understanding the transient X-ray sources. 

In all the cases, the ion and electron temperatures merge or tend to merge at around transition radius. This is because, the electrons are 
in thermal equilibrium with the ions and thus virial around the transition radius, particularly when Mil. As the sub-Keplerian flow advances, 
the ions become hotter and the corresponding temperature increases, rendering the ion-electron Coulomb collisions weaker. The electrons, on 
the other hand, cool down via processes like bremsstrahlung, synchrotron emissions etc. keeping the electron temperature roughly constant 
upto very inner disc. This reveals the two temperature flow strictly. 

Important point to note is that we have assumed throughout the coupling between the ions and electrons is due to the Coulomb scattering. 
However, the inclusion of possible nonthermal processes of transferring energy from the ions and electrons (Phinney 1981, Begelman & 
Chiueh 1988) might modify the results. However, as argued by Narayan & Yi (1995), the collective mechanism discussed by Begelman & 
Chiueh (1988) may dominate over the Coulomb coupling at either a very low a or a very low M. Instead, the viscous heating rate of ions 
is much larger than the collective rate of nonthermal heating of electrons, unless a is too small what we have not considered in the present 
cases. Therefore, the assumption to neglect nonthermal heating of electrons is justified. 

Now the future jobs should be to understand the radiation emitted by the flows discussed here and to model the corresponding spectra. 
This will be the ultimate test of the model in order to explain observed data. 



APPENDIX A: DISCUSSION OF BOUNDARY VALUES 

We have four coupled nonlinear differential equations (6), (10), (12), (17) to be solved for &, c s , A, T e ; equations also involve p and P. 
To eliminate p and P, we use mass transfer equation (1) and equation of state (9). Therefore, in total we have five differential equations 
supplemented by an equation of state. Hence, we need five boundary conditions to start integration. Equation (1) can be integrated to obtain 
M already given in Eqn. (2), which is supplied as an input parameter. Similarly, integrating Eqn. (10) we can obtain angular momentum flux 

M(X-X^) = -AnAW^\, (Al) 

where A m is the specific angular momentum at the inner edge of the disc, to be fixed by no torque inner boundary condition. Note that A m < A c 
(see, e.g., Chakrabarti 1996). 

We therefore need the initial values of &, c s , A and T e to solve the set of equations. When we impose the condition that the flow must 
pass through a critical radius x c (around a sonic radius) where D = 0, & and c s at x c are related by a quadratic equation of Mach number 
given by Eqn. (24). 

For the continuity of d&/dx, N = at x c . Therefore, from N = which an algebraic nonlinear equation, c s at x c can be computed 
iteratively (using bisection method), which in turn fixes & at x c as well, provided A is known at that radius. 

Now we need to set appropriate values of x c and corresponding specific angular momentum A c . This is fixed iteratively by invoking the 
condition that the critical point to be saddle-type. This can be seen as follows. First we impose that A m = A c . Then by fixing the value of A c we 
find that if x c is greater than a certain critical value x cc , then the type of x c changes from saddle-type to 0-type which matter never can pass 
through. Figure lATT i shows how the type of critical point and then solution topology change with a slight increase of x c . On the other hand, 
as x c decreases from x cc which corresponds to the energy at x c (E c ) increases, the sub-Keplerian disc decreases in size. This advances the 
Keplerian disc. The reason is that increasing E c corresponds to decreasing x c and then increasing centrifugal energy (A 2 /2x\ c ) which keeps 
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Figure Al. Comparison of the variation of radial velocity as a function of radial coordinate (a) between solutions with x c = 5.5 (dotted curve) and x c = 5.7 
(solid curve), when A c = 3.2, (b) between solutions with A c = 3.2 (dotted curve), x c = 3.3 (solid curve), when x c = 5.5. Other parameters are same as that in 
Fig.|2]forM = 0.01. 



the flow Keplerian until inner region. However, in principle the solution of the model equations is possible to obtain for any value of x c from 
x cc to the marginally bound orbit xt. 

Once x c is fixed at x cc , we have to obtain the best value of A c . By increasing the value of A c beyond a certain critical value A cc at a 
particular x c , we again find a transition from saddle-type to O-type critical point. Fieure lATb shows how the type of critical point and then 
solution topology change with a slight increase of A c . On the other hand, decreasing A c from A cc will tend the disc to more Bondi-type. Now 
for A cc we again have to obtain a new value of x cc following the procedure outlined above and thereafter corresponding A cc . This needs to be 
continued iteratively until a specific combination of critical radius x cc and corresponding specific angular momentum A cc , lying in a narrow 
range, is obtained which leads to a physically interesting large sub- Keplerian accretion disc, when matter infalling from a largest possible 
transition radius to a black hole event horizon through a saddle-type critical point. However, in principle the solution of the model equations 
is possible to obtain for a range of A c such that A cc > A c > 0. 

There is, however, another initial value namely T e at x c (T ec ) to be assigned. Choice of T ec depends on the observed nonthermal radiation 
which restricts the value of T e in general. But this restriction can only provide an order of magnitude of T e . An exact value of T ec should be 
obtained iteratively from a plausible range of T e at x c so that the values of x cc and A cc obtained following the above mentioned procedure 
converge. 



APPENDIX B: COMPUTATION OF DERIVATIVE OF VELOCITY AT THE CRITICAL RADIUS 

We first recall the derivative of velocity from Eqn. (20) 

M = N(x,&,c s ,A,T e ) 
dx D(&,c s ) ' 

At the critical radius 
d& 
dx 

Therefore, we apply l'Hospital's rule and obtain 

d& £[N(x,&,c s ,A,T„)] 
dx ~ &VW,c.)1 

dN_ , ,JN_M + dN_dij_ + M dA | dN AT C 

dx d& dx dc s dx d.X dx dT e dx tTt1\ 

~ dD d& ! dD dc s ' ' 

dft dx dc x dx 
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Now combining with Eqns. (26), (27), (28) we obtain 



d§ Ni + N 2 f 



dx dc, c, dX 



such that 

dd 
dx 



N 2 -D 1 ± V(»i - N 2 ) 2 + 4D 2 Ni 



2D 2 

where upper and lower signs correspond to wind and accretion respectively. 



(B4) 



dx D[ + D 2 f 

where 

dN dN J dN jc 2 -2xaj \ dN (T, - 1)4tt c,x 3 ' 2 , . \ (J 



dN dNjc s &\ dN I lax I„ +l I c] \ \ dN & 

N ^M + dF s (l--) + ll(^-[i- + a j + W} 1 ~ r ' )r ^' (B6) 

dD J 

», = —- (B7) 

dc s c s 

dD dDlc. &\ 

D 2 = — + — — ) . (B8) 



dft dc s \ i? c s 
Finally cross-multiplying in Eqn. <B4t we obtain a quadratic equation 
(d&\ 2 (d&\ 

D2 Ux] + (Dl ~ m W~ Ni = (B9) 



(BIO) 
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